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: ABSTRACT 

Context. Comet 67P/Churyumov-Gerasimenko is the main target of ESA's Rosetta mission and will be encountered in May 2014. As 
^SJ ' the spacecraft shall be in orbit the comet nucleus before and after release of the lander Philae, it is necessary necessary to know the 

conditions in the coma 

Aims. Study the dust environment, including the dust production rate and its variations along its preperihelion orbit. 
Qh ' Methods. The comet was observed during its approach to the Sun on four epochs between early- June 2008 and mid- January 2009, 

' over a large range of heliocentric distances that will be covered by the mission in 2014. 

Results. An anomalous enhancement of the coma dust density was measured towards the comet nucleus. The scalelength of this 
enhancement increased with decreasing heliocentric distance of the comet. This is interpreted as a result of an unusually slow expan- 
sion of the dust coma. Assuming a spherical symmetric coma, the average amount of dust as well as its ejection velocity have been 



^ ' derived. The latter increases exponentially with decreasing heliocentric distance (rn), ranging from about 1 m/s at 3 AU to about 25-35 

I m/s at 1.4 AU. Based on these results we describe the dust environment at those nucleocentric distances at which the spacecraft will 

, presumably be in orbit. 



Key words, comets: general — comets: individual (67P/Churyumov-Gerasimenko) 



' 1 . Introduction are very faint when far from the Sun, and require observations 

. „„., ^ ,. , , with large telescopes for which observing time is difficult to re- 

In May 2014 ESA s Rosetta spacecraft, an orbiter and Ian- ^^.^^ r^^.^ j^^^ particularly to the short-period comets of the 

rn der mission to a comet, wil encounter its target, comet jupjter-family, of which the Ro setta target com et 67? is a mem- 

O 67P/Churyumov-Gerasimenko hereafter 67P), at a heliocen- j ^ ^ -^^^^ (lA'Hearn et aL L [1981 is a ques- 

, tnc distance rj, of about 4 AU. It will go in orbit around ^.^^^^^^ ^^^^^^^ ^^^^ production rate because the condi- 

O ■ the nucleus in September 2014, when the comet is at rh= 3.4 ^.^^^ ^.^^^^jj ^^^^^^^^ ^ ^ ^.^^j^^ (homogeneous, isotropic 

^ ; Mi and start global observations and mapping of the nucleus ^^^^^^^^ ^^^^^^ ^j^^ '^^^^^ f^j^jj^^ 

T-H . jGlassmeier et al., 2007) . At rj,- 3 AU, the lander Philae will ^^^^^ ^^^^ expansion velocities of the dust in the coma are 

^ land on the surface and perfom the first-ever in-situ analysis of j^^^ ^^^^j^^^^ ^ ^^^^j ^^^^ ^^^^j^^ ^^^^^j^^^ 

.r^ . comet nucleus material. The orbiter will be monitoring the evo- .^^ ^^^^^^^ ^^^^^ paiticulai'ly for the larger grains. 

. lution of the nucleus and the coma along the comet s pre- and ^^.^ j^^^ ^ . j^, .^.^^^ ^^^^^ reflecting the 

H post-perihehon orbit for more than one year while approaching ^^^^^ just production of the nucleus over a range of solar 
. . . the nucleus as close as a few km above the surface. 



When approaching the Sun, at distances from about 5 to 2 
AU, cometary nuclei usually switch on the activity and produce 
a gas and dust coma characterized by a changing environment 
due to secular and short-term effects, e.g., solar heating and ro- 



distances. 



Apart from the scientific interest in obtaining a better knowl- 
edge of the dust coma environment and evolution of 67P along 

„ , , ^ , • J- , J the preperihelion orbit, this knowledge is required for the safety 

tation of the nucleus. Consequently, an increase of the produc- . , , • r , r, ■ ■ *^ , 

, . J r , J ■ J 1 and the operations planning of the Rosetta mission. After the 

tion rate and expansion speed of the dust is expected as a result „ i u j i • ^rv/?o ^ir, i . j 

„ , . J 1 .c 1 1- • • / 1 1 KOiefta launch delay in 2003, 67P was selected as the new tar- 

of the increased heating of sublimating ices from the nucleus 



Unfortunately not much is published on the coma status and the 
dust environment of comets at large distances, because comets 



get of the mission replacing the original target 46P/Wirtanen 
which was much better characterized by Earth-based observa- 
tions. At the time of its selection 67P had just passed its perihe- 
* Based on observations collected at the European Organization for lio"- Thus, most observations existing to date were obtained dur- 
Astronomical Research in the Southern Hemisphere, Chile (programs ing t he subsequent post-perihelion phase (Schulz et al. , 2004|; 



381.C-0123 & 082.C-0740) and Telescopio Nazionale Galileo (TNG) Wei ler et al. L 120041) . wh ereas only a few observations exist of 
of INAF (program TAC.35) the pre-perihehon phase (ISchleicher I l2006l iLamv et all l2006l 
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Table 1. Summary log for the observations of comet 67P. The 
table lists the date of observations (YYMMDD), the heliocen- 
tric rh and geocentric distance A of the comet in AU, the phase 
angle (a), and the position angle of the extended Sun-Comet ra- 
dius vector (PA,fl,7) and of the velocity vector (PAvei.)- The last 
column gives the multiplicative facto r to apply to transf orm the 
measurements to a phase equal to 0° d Schleicher L |20 1 Ol) 



Date 


r;, 


A 


a 


PA„,7 




CoefF. 


YYMMDD 


AU 


AU 










080601 


2.98 


2.50 


18.83 


252.9 


250.3 


1.89 


080604 


2.96 


2.44 


18.67 


253.3 


250.2 


1.89 


080605 


2.95 


2.42 


18.60 


253.4 


250.2 


1.88 


080904 


2.30 


1.40 


14.30 


55.7 


260.2 


1.67 


080906 


2.28 


1.39 


15.23 


57.4 


260.3 


1.71 


080908 


2.27 


1.39 


16.14 


59.0 


260.6 


1.76 


081022 


1.93 


1.51 


30.63 


73.0 


259.0 


2.44 


081026 


1.90 


1.53 


31.37 


73.3 


258.1 


2.47 


090113 


1.36 


1.67 


35.98 


68.7 


239.1 


2.64 



I2OO7II2OO8I) . except from thos e obtained by amateu r astronomers 
(cara.uai.it, www.aerith.net), ('K idgerLl2003ll2004 . 

The 2009 return of the comet was the first and last oppor- 
tunity to measure the pre-perihelion dust environment of 67P 
before Rosetta will have arrived at the comet. 

Here, we report the results of our dust coma analysis of 
67? using observations obtained during the pre-perihelion phase 
2008-2009 between rh~ 3 to 1 .4 AU. Our obse rvation s are com- 
plementary to those reported bv Tubiana et al. I (*2008. 20 10|), ob- 
tained at larger heliocentric distances and those by Lara et al. I 
(I2OIOI) . obtained during the post-perihelion phase. 

2. Observations 

Comet 67P was observed at four epochs (see Table [1} dur- 
ing the inbound phase before it reached its perihelion in Feb 
2009. During the first three epochs the observations were ob- 
tained at the 8.2m Very Large Telescope (VLT) of the European 
Southern Observatory in Chile using the FORS2 focal reducer 
instrument (see www.eso.org/sci/facilities/paranal/instruments). 
The last run, when the comet was in the northern hemisphere, 
was performed with the Telescopio Nazionale Galileo TNG 
at La Palma using the camera-spectrograph DOLORES (see 
www.tng.iac.es/instruments/). During each run a series of im- 
ages was acquired using broadband filters {VRI bands) and 
sometimes narrower band filters, taken a few days apart in order 
to check the comet short term variability. Differential tracking 
of the telescope at the cometary velocity was applied for all 67? 
images. Since the comet was in front of star-rich regions and 
since it had a relatively fast proper motion, each run consisted 
of several (5-9) exposures in order to reduce the star contribu- 
tion by median averaging of the images during data processing. 
At the VLT the observations were performed in service mode 
and the flux calibration relies on the zero point estimations (Zp) 
provided by the observatory for each photometric night. The ob- 
servations at the TNG were taken in visitor mode and a series of 
photometric standards were obtained for calibration purposes. 

3. Data reduction 

All images were processed using standard CCD reduction proce- 
dures (BIAS & Flat Field correction). Then a first-order constant 



sky level, measured in image regions away from the comet's 
photocenter, was subtracted. After combining all filter frames 
of a single night (aligned to the comet and median averaged), 
possible sky residuals were checked in the resulting frame and 
subtracted using the EA/ function (see below). Finally, the im- 
ages of each night were calibrated in A f, i.e. t he albedo multi- 

flied by the filling factor of the dust in the coma (lA'Hearn et al. L 
984; Tozzi et al. , 2007). In order to allow comparison of the 
results obtained at different phase angles, the calibrated images 
were re-calibrated for the phase equal to 0° using the multiplica- 
tive phase correctio n coefficient (last column of Table [1]) given 
by Schle icher et al I (11998) and recently refined in ( Schleicher L 
2010). The observation series taken over a few days apart were 
searched for the presence of possible short-term variability. 
Since day-to-day variability was not found, all images obtained 
with the same filter during each epoch were median averaged, 
thus further reducing the contribution of the background stars 
and increasing the signal-to-noise ratio of the cometary coma 
flux. 

These final images in the VRI filters allowed to assess the 
dust distribution in the coma of 67P. The gas contamination in 
the broad band filters was negligible for all four epochs, as was 
verified in spectra and measurements through an interference fil- 
ter with the passband centered in regions without any gas emis- 
sion (central A - 8340 A, AA - 480 A). No or neghgible gas 
emission was present in this relatively wide passband, which was 
verified a posteriori also by the spectra. 

4. Data analysis 

The final resulting images were initially analyzed using the 
2A/(p) function. lA f(p) is proportional to the average column 
density of the solid component at the projected nucleocentric 
distance p. It is equal to InpAfip) and is measured in cm. As 
shown in i Tozzi et alTI (2001), l,Af(p) is constant with the pro- 
jected nuclear distance p in the case of a comet with a dust out- 
flow of constant velocity and production rate, and if sublimation 
or fragmentation of the grains are excluded. Only the solar radi- 
ation pressure introduces a small linear dependence, but usually 
only at large distances from the nucleus, definitely larger than 
our field of view (FOV). Unfortunately, despite of all efforts, our 
final images suffered from remnant flux contamination of the nu- 
merous background stars so much that the stars were not com- 
pletely erased even after using the median average. However, it 
was always possible to check and verify the radial profile of the 
SA/ function up to about 10^ km from the nucleus. Figure [T] left 
panel, gives an example of a measured ZA/ profile. Clearly, it is 
not constant with radial distance p from the nucleus and shows 
a fast decrease as p increases, until it reaches a constant value 
at about 20000 km. The A fp function, shown on the right panel 
of the same figure, is also strongly dependent on p and clearly 
cannot be used as a proxy of the dust production rate. 

Using a trial & error procedure the profiles were reproduced 
by integrating along the line of sight of a spherical symmet- 
ric coma with an "optical density distribution" of the dust as 

N{d) = [Nc + Me"*^* + N^e'^i'']l(f- where d is the (non pro- 
jected) nucleocentric distance in km. This function has a con- 
stant term A^c representing the constant optical density and two 
exponential terms with A^s and A^i describing a short-scalelength 
and a long-scalelength optical density increases, respectively. 
The cf' term represents the spatial density attenuation due to ge- 
ometric expansion of the dust in the coma, L\ and Ls are the 
relative scalelengths of the long and short-scalelength variations 
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p (km) p (km) 

Fig. 1. Left: Measured R filter EA/ in function of the projected nucleocentric distance p for observations at rh~ 2.3 AU together 
with the model functions. Crosses represent the measurements, the solid line shows the total model function, the dashed line gives 
the part related to Nc+Nj and the dotted line that of A^s- The "bumps" (for instance, that at about 32000 km) are due to the remnant 
flux from background stars, easily identifiable by eye in the images. Note that the profile has been modeled up to 10^ km, but here 
a smaller part is shown. Right: Measured Afp function for the same date and filter This shows that the Afp is strongly dependent 
on p and consequently it cannot be used as proxy of the dust production 



in A^(d). Nc, N\, Ns describe the "optical density distribution", i.e. 
the albedo multiplied by the light scattering area of the grains in 
a cubic centimeter They are measured in cm"' (area over vol- 
ume, cm^/cm^). Of course, in order to derive the corresponding 
spatial density of the dust (number of grains per cubic centime- 
ter) one has to consider the albedo and size distribution of the 
grains. 

For comparison with the data, A^(d) is integrated along the 
line of sight. Figure [T]provides an example of the computed pro- 
file as a function of the projected nucleocentric distance (p) us- 
ing the adopted model function A^(d). The agreement between 
the measured radial profiles and the models, defined above, is 
very good. The only differences are within the "bumps" caused 
by the background stars. The best values of the model parameters 
A^c, N\, Nfi, and L\ and Ls together with the observing geometry 
and filters are listed in Table |2]for the four observing epochs. 
The errors indicated in the table refer mainly to the dominat- 
ing systematic errors. is the most affected by the errors due 
to contaminating flux from background stars. Our checks of the 
lAf profiles up to more than 10^ km, identifying the star con- 
tribution visually, indicate an uncertainty of Nc not greater than 
a factor two. A^s estimation is also affected by the systematic er- 
rors because is of the order of the seeing, in particular for the 
observations taken during the first epoch. Hence, for that epoch, 
Ls should just be considered as an upper limit and, consequently, 
A^s as a lower limit. The parameters less affected by the system- 
atic errors are N\ and L\, because L\ is always much larger than 
the seeing and the intensity of this component is high enough so 
that the contribution of possible background stars is negligible. 
By varying these parameters in the adopted function for A^(d), 
we estimate that the uncertainties are better or of the order of 
25% for M andLi. 

The last three columns of Table |2] are the optical cross 
sections (SA), measured as the integral of A^(d), of three 
components: constant(SAc), long-scalelength (SA;) and short- 
scalelength(SAs). They represent the total area covered by the 
particles of each component multiplied by the particle albedo. 
At first order, SA is independent of the seeing and gives "real" 



value (not just upper or lower limits) also for the short scale- 
length component. For the constant component, we arbitrarily 
computed its optical cross section assuming a limit for the nu- 
clear distance of 5 x 10^ km, since in principle it extends to infin- 
ity. The results show that, while SA of the constant component 
does not change much with rh, the SA values of the other two 
components change a lot (by factors up to 20-40). We conclude 
that the two exponential components were replenished by fresh 
grains produced as the comet approached the Sun. However, we 
cannot exclude a small replenishment of fresh grains also in the 
constant component, because the grain ejection velocity should 
increase with decreasing rh. 

It is important to note that the flux contribution from the nu- 
cleus was not negligible when the comet was at rh~ 3AU. With 
the radius = 2.0 km and the geometric albedo i n the R filter A 
= 0.054 (iLamv et al. L l2006l 120071: iKellev et al.L l2009h . the op- 
tical cross section (SA) of the nucleus at 0° phase angle is equal 
to 0.68 km^. Taking in to account the facto r 4 due to the different 
definition of albedo (lA'Hearn et al. L[T984l) . the optical cross sec- 
tion of the nucleus corresponds to 2.7 km^, which is a large part 
of the measured SA, = 3.9 km^ of the short-scale length cross 
section in our R filter For the other epochs, the contribution of 
the nucleus is always negligible. 

In the past, such anomalous enhancements were observed in 
other comets and interpreted by the presence of organic grains 
that sublimated in the coma with a certain lifetime t w hile mov- 
ing away from the nucleus (iTozzi et al. L l2004 l2007h . If this is 
the case the lifetime (and the scalelength) should decrease when 
the comet approaches the Sun. Our results for 67P show an in- 
crease of L\ with the comet approaching the Sun, which ex- 
cludes that grain sublimation or any other photolytic process 
were responsible for the enhancement. 

Figure |2] gives the scalelength Li as a function of the he- 
liocentric distance rh. The scalelength values are well fitted 
by the function function L\= crj^^ with best fit given for c = 
(195000 + 29000) km and yi ='(-3.93 ± 0.20) in V and c = 
(120000+ 18000) km and n = (-3.20 + 0.19) in 7? . 
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Table 2. Fit parameters of the lAf profiles. The first two columns give the average date of observations and the filter used (V , R , 
I ). Nc, Ni and A^s are the optical density distributions of the constant, long and short components, respectively. Li and Lj are the 
respective equivalent scalelengths. The last three columns give the total optical cross sections of the three components. 



Date 


Filter 


N, 


A'l 


u 








SA, 


SA, 


SA, 


YYMM 






lO^'^cm-' 


km 




km 


km^ 


km~ 


km^ 


0806 


V 


02+" "^ 


0.30+0.07 


2930+ 


60 


0.42±0.10 


793±16 


14±^4 


8.0+2.0 


3.0+0.6 


0806 


R 


02+'' ''^ 

"•"^-0.01 


0.19±0.05 


3750+ 


75 


0.40±0.10 


1056±21 


14±^'' 


6.5±1.6 


3.9±1.0 


0806 


I 


02+°-°- 
»■!!! 


0.31±0.08 


3170+ 


63 


0.25±0.06 


800±16 


15±''' 
30±° 
24±|^ 


9.0+2.2 


1.7±0.4 


0809 


V 


04+'' '''* 
8-Q3 


0.31+0.08 


7350+ 


145 


0.99±0.25 


976+20 


20±5 


8.9±2.2 


0809 


R 


03+™ 

"•"-'-0.015 


0.26±0.06 


9300+ 


186 


0.74±0.18 


1330±27 


23±6 


9.3±1.8 


0809 


I 


03+™ 

"•"-'-0.015 


0.26±0.06 


9600+ 


190 


0.77±0.15 


1390±28 


23+23 

^-'-12 


23±6 


10.3±2.6 


0810 


V 


06+""^ 

"•""-0.03 


0.26+0.06 


13000+ 


260 


2.08±0.52 


1300+26 


44±22 


32±8 


26±6 


0810 


R 


13+°" 

"•^-'-0.06 


0.53±0.13 


13000+ 


260 


3.43±0.61 


1380±28 


92+'2 
79^^^ 

"-38 


66±16 


45±11 


0810 


1 


11+"" 
"•^^-0.05 


0.50±0.12 


11300+ 


226 


2.82+0.70 


1380±28 


54±13 


37+9 


0901 


V 


"•"-"-0.02 

03+™ 

"•"■^=^0015 


0.42±0.10 


62000+ 


1240 


2.38±0.60 


4300±86 


40+^*" 


282±70 


96±24 


0903 


R 


0.41+0.10 


46600+ 


930 


2.51±0.63 


4700±94 


206±51 


111±28 



E 




<u (AU) (AU) 



Fig. 2. Measured L\ for V (star) and R (square) filters as a func- 
tion of the heliocentric distance. The lines show the best fit as- 
suming the function described in the text. Solid line represents 
the V and dotted line the R filter. 



Fig. 3. Computed grain ejection velocity as function of the 
heliocentric distance for observations in V (solid line) and in 
R (dotted Une). 



If we assume that this long-scalelength component was due 
to fresh dust produced in a spherically symmetric coma during 
the approach of the comet to the Sun, the coma was expand- 
ing with a velocity given by which is equal to c x yi x 

('"/i(O)'-^'^'' X Since in this part of the orbit the hehocen- 
tric distance varies linearly with time as rh = r/,o + voTp with Tp 
- days to perihelion, r/,o = 1.011 AU and vq - -7.18 x 10"^ 
AU/day, the time variation of the scalelength is ^^^p- = c x vq x 
7i X (r/,(f))^^'"''. This is the expansion velocity of the grains 
that are at nucleocentric distance equal to L\. It is known that 
the grains are accelerated by gas drag only in the inner part 
of t he coma, i.e. within few radii of the comet nucleus (see 
e.g. ICombi et al. I (Il997h ^ and that they move with almost con- 
stant velocity beyond that limit. They are affected only by solar 
radiation pressure, efficient at large scales, and by the gravity of 
the comet, efficient only very close to the nucleus. This means 
that, roughly speaking, the grains expand in the coma with a 
widely distributed, but unchanging velocity. Hence, the grains at 
nucleocentric distance L\ were ejected some time before (AT") 
with the velocity To compute Ar we have just to divide Li 
by ^^^It^. This gives AT" - -^^^ x rh(f), i-e. the time necessary 



to reach L\ is linearly dependent on rj,. For example, for obser- 
vations in the R fflter, Ar is about 87 days at rh= 2 AU. The 
value of ^^p- for the same rh and fflter is about 150 km/day or 
1.73 m/s, i.e. the expansion velocity is quite low. As the grains 
are expanding almost freely (see above), the expansion veloc- 
ity at a certain time Ti is equal to the ejection velocity at time 
Ti - AT. By solving the equations numerically, we get the ejec- 
tion velocities as functions of the heliocentric distance, as shown 
in Figure|3] They are equal to Vej = Ve/o x rj^^ ' = Vejo x rj^, with 
Vejo - 195 ± 21 and 100 ± 1 1 m/s in V and R , respectively. 

This formula is similar to the ejection velocity assumed in 
models, i.e. Vej ■\/^ x rj^, with /3 = the ratio of the radi- 
ation pressure over solar gravity forces (see e.g. iFulle et aTl 
(l2010h V We get y = -(4.93±0.20) and -(4.20±0.20) in V and 
R , respectively, which is quite different from the assumed 
values: -0.5 (Ishiauro , 2008; Kelley et al. , 2008, 2009) or -3 
(Agarwa l et al. 1120 10). The ratio of /?, of the particles observed 
in V with respect to those observed in R is about 3.8± 1 .2, a rather 
high value that is difficult to explain. 

It is important to note that all results about the grain ejec- 
tion velocity assume isotropic outffow of the dust. If the comet 
showed an asymmetric emission (e.g., a jet) in the direction of 
the observer, this would have changed the results leading to a 
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higher ejection velocity. However, this would mean that the jet 
direction and intensity had to remain unchanged for about four 
days during each epochs, since even the images taken several 
days apart during each epoch did not show any change in flux. 
Such a stable jet could only be produced if the spin axis of the 
nucleus was pointing exactly into the observer's direction, which 
could be the case only for one of the three epochs if at all, but 
certainly not for all three epochs. 

5. Discussion 

5.1. Other Observations 

It is important to note that during the post-perihelion phase 
comet 67P displayed a completely different behavior from what 
we found in the pre-perihelion images. The lAf profiles did 
not show any anomalous enhancement towards the nucleus. 
For example , the pr ofiles obtained from HST observations (see 
iLamv et al~l ( l2006l) ) and downloaded from the HST archive, 
show only the signature of the nucleus itself at p close to zero, 
but are otherwise 'flat' up to p = 4 x 10'* km, the FOV of the im- 
ages. A similar behavior is also found in other post-perihelion 
measurements, which demonstrates that the slope parameters 
of log(I) ys log(p ) are equal to about -1 (iLara et al. 
I S chleicherT l2006h . Obviously, after perihelion the comet had al- 
ready developed an extended coma with the conditions necessary 
for a reliable determination of Afp (see above). 

5.2. Grain size 

As shown above, for heliocentric distances greater than ^ 2 
AU, the dust expansion velocities in the coma of 67P were 
about two orders of magnitude smaller tha n those usually 
claimed for cometary c omae (see Crifo and R odionovl (Il997h . 
[Foster and Green. 1 (1200 7) and references therein ). This may be 
an indication that the grains in 67P were big. lAgarwal et al. I 
(l2007l 1201 Ol) computed that velocities of the order of 30-50 
m/s at rh= 1.3 AU correspond to grain radii ranging from about 
10 yum to 1 mm (depending on the emission scenario) assum- 
ing spherical grains with a density equal to 1 g/cm^. Since it is 
likely that the grains are fluffy and non spherical (and the density 
much lower than than a ssumed), their size can be much larger. 
iHadamcik et al. I (l2010l) have also reported an anomalous slope 
log(I) vs log(p) (close to - 1 .5) in the inner coma (p between 2000 
and 8500 km) in December 2008, which they interpreted as be- 
ing due to the presence of large particles. The presence of large 
and fluffy particles is also consiste nt with polarization me asure- 
ments performed in March 2009. IHadamcik et al. I (12010 ) note 
that the high polarization ( 6%) observed near the nucleus is 
ty pical for small, sub-micron particl e s (cf. comet Hale-Bopp 
in IHadamcik and Levasseur-RegourdI (l2003h '). This contradic- 
tion can be resolved by suggesting large, but por ous particles 
(Hadamcik et al. ','2010*). Note that several authors (Fulle et al. , 
2004, 2010; Moreno_elaL_, ,2004; Kelley et al , 2008; Ishiguro, 
2008j)) suggested the presence of grains in the cm-range based 
on the observations of a trail, the tail and necklines in 67P. 

5.3. Grains Density in fhe Coma 

For regions with <«; <K Li, which are relevant for the 
Rosetta spacecraft during its close approach or when orbiting 
the nucleus of 67P, the mean optical density distribution of the 
dust at a nucleocentric distance d is simply equal to A^(d) - 
(Nc +Ni+Nf,)/d^ - N{)ld^, with No represents the optical density 



distribution std - I km. For example, at 1 .93 AU from the Sun, 
we find M) = 4.1 x 10"** cm"' (see the R filter results in Table|2). 
Assuming the power-law size distribution n{r) - iidij-y'^, equal 
for the three components of A^(d), the optical density is equal to 

with A = the grain albedo, and the integral extended from the 
minimum (r,„) to the maximum (r^) radius of the grains as- 
sumed to be spherical. The local density mass due to the dust is 

with cr being the average density of the grains. The ratio is 

then independent on d and is equal to | j y3-t_^-t-^ ■ Assuming 
Tv '"m and A: < 3, as big grains dominate, we have M(d) = 
A^(d) I J l^f Ti - So, for the given physical parameters of the grains 

(albedo and density), the ratio depends on the size of the 
biggest grains in the dust size distribution, r^ and, weakly, on 
k. With a grain density cr = 0.2 g/cm^, an albedo A - 0.04, 
and k - 2 the mass density is M(d) = 3.3N{d)r^. For exam- 
ple, for rh~ 2 AU, at 10 km from the nucleus, the optical den- 
sity is 4.1 X 10"'"cm"'. Assuming a maximum grain size of 1 
mm, the mass density would be 1.4 x Hence, dur- 

ing one orbit revolution at this distance each square meter of the 
Rosetta spacecraft would intercept a total mass of only 8.6 g. 
If the maximum grain radius is 5 cm the mass density would be 
6.5 X 10"'g/cm^ and the mass intercepted would be 430 g, a non- 
negligible amount, considering that the maximum cross section 
of the spacecraft is 17 m^. 

It is important to note that the main contribution to the den- 
sity of the coma in regions close to the nucleus is produced by the 
short scalelength component, that sometines is affected by the 
seeing of the observations. In this case N j may be a lower limit. 
However, for the considered example at rh~ 2 AU this compo- 
nent is not strongly affected, because its scalelength corresponds 
to about 1.3", which is larger than the seeing of the observational 
nights. 

6. Conclusions 

The anomalous enhancement of the dust density of the coma of 
67P towards the nucleus, measured in 2008 and 2009, can be 
explained by a very slowly expanding dust cloud, with veloci- 
ties of the order of 1 m/s at rh=3 AU. The slow dust expansion 
velocity supports a scenario where at large pre-perihelion helio- 
centric distance the coma of 67P was dominated by large grains, 
with dimensions greater than 10 jum - 1 mm, if they are spherical 
with a density equal to 1 g/cnr', or much larger, if they are fluffy 
aggregrates. By modeling the radial coma flux profiles with a 
spherically symmetric coma, the optical density distribution of 
the dust was quantified as a function of the nucleocentric dis- 
tance d. The ejection velocity was derived as a function of the 
heliocentiic distance, rh. The model also allowed to derive the 
mass density of the dust as a function of d assuming a grain size 
distribution in which big grains dominate the coma. 
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